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Introduction 

Quantum Chromodynamics (QCD) at finite tempera- 
ture and density is a very active area of research. The 
equation of state of QCD and the nature of the transition 
from the hadronic phase with broken chiral symmetry to 
the chirally symmetric, deconfined quark-gluon plasma 
phase is of great importance for a better understanding 
of the experimental data, e.g. [1]. 

For full QCD with dynamical quarks one expects 
the deconfinement phase transition to be a crossover as 
quarks expUcitly break the underlying center symmetry 
of the gauge group. The nature of the chiral phase transi- 
tion primarily depends on the value of the current quark 
mass, which explicitly breaks chiral syinmetry, as well as 
on the strength of the chiral anomaly [2]. While the con- 
finement phase transition is driven by gluodynamics, the 
chiral phase transition is governed by strong quark inter- 
actions. It is a highly non-trivial observation that both lie 
remarkably close at least for small quark chemical po- 
tentials [3, 4]. An understanding of this interrelation is 
subject of an ongoing debate. 

Hence, the resolution of the QCD phase diagram re- 
quires both, the computation of the relevant observables, 
in particular the order parameters, and also some analytic 
understanding of the mechanisms involved, see e.g. [5]. 
Such a twofold task is best obtained within a combina- 
tion of different methods that allow a direct access to the 
observables as much as to the mechanisms behind the ob- 
served phenomena. Functional continuum methods such 
as functional renormalisation group equations (FRG) and 
Dyson-Schwinger equations (DSE) are well-suited for 
the above task. They are indeed complementary to lattice 
simulations: most importantly they are amiable towards 
chiral fermions including the anomalous breaking of chi- 
ral symmetry and are straightforwardly appUed to finite 
chemical potential. Ideally lattice simulations and func- 



tional approaches go hand in hand, and should be used 
to improve and to check respective results, also leading 
to better systematic error estimates. In combination, this 
should allow us to map out the phase diagram of QCD. 

Functional RG 

We shall present results on the order parameter for the 
chiral and the confinement-deconfinement phase transi- 
tion obtained with functional methods, mostly with the 
FRG. Within a functional approach one usually com- 
putes correlations functions of quarks and gluons or that 
of composite operators such as hadronic degrees of free- 
dom. These correlations functions are related by an in- 
finite hierarchy of partial integro-differential equations 
which are solved within specific approximations to the 
full system at hand. Clearly a sound discussion of the ap- 
proximations is of chief importance for the reliability of 
functional approaches. However, it is beyond the scope 
of the present overview, and for more details we defer 
the reader to the original works. 

Most of the results presented here are obtained in 
Landau gauge QCD with the classical action 

5qcd = \jF^^,'' + ^j{dA-f + jc%dDr'& 

+ j ^r{iD + im^ + iHYo)Vf, W 

where — > 0. An infrared regularisation is achieved with 
the introduction of momentum-dependent mass-terms 
for ghost, gluon and quark fields. We also introduce cut- 
off terms for effective low energy degrees of freedom 
such as mesons and baryons. This leads to a scale- 
dependent effective action, r;t[A,C,C, VA,i//;(7,7r,...], 
with infrared cut-off scale k. The a and n fields stand for 
mesonic composite operators. For A: — )• 0, we approach 



the fully interacting theory, whereas for A; — ^ oo we are 
left with the asymptotically free high energy QCD. 
An infinitesimal change of the scale k is described by 
Wetterich's flow equation, [6], 



kd. 



d^nm = ^Tr 



(2) 
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where t =lnk/A, and = {A,C,C,\jf,\jf;a,7i,.. 

(2) 

nally, FJ. denotes the second derivative of F^. w.r.t. the 
fields. The cut-off functions Rj^^ provide infrared cut-offs 
for all fields, including the composite fields, i.e. a and n, 
for details see the reviews on gauge theories, [7, 8, 9, 10]. 
We emphasise that (2) comprises a first principle QCD 
flow, the appearance of composite operators does not sig- 
nal an effective field theory setting but rather a conve- 
nient parametrisation, see [8, 11]. The flow (2) has a sim- 
ple diagrammatic form depicted in Fig. 1, see [12]. The 



FIGURE 1. Functional flow for QCD: the lines denote full 
field dependent propagators. Crosses denote the cut-off inser- 
tion d,R. 

first and second loop generate gluon and ghost fluctua- 
tions respectively, the third loop generates fluctuations 
of the quarks, and the last loop stands for the loops of 
the mesonic a- and tt -fluctuations, and possible further 
hadronic degrees of freedom. 

Flows for Yang-Mills propagators in the Landau gauge 

The Landau gauge has very peculiar properties that fa- 
cilitate the computation of correlations functions: first of 
all, the ghost-gluon vertex is protected from renormal- 
isation. Second, we have infrared ghost dominance. Its 
weak form relevant for the present investigations sim- 
ply entails that the gluon dressing function is vanishing 
in the infrared, p^/r^^\p — > 0) — > 0, whereas the ghost 

dressing function does not, p^ /Y^^\p — 0) > 0. For a 
detailed discussion of the global properties of Landau 
gauge Yang-Mills theory see [13, 14]. A related gauge 
with similar properties is Coulomb gauge, for a FRG 
study and further references see [15, 16]. 

The full set of flow equations for ghost and gluon 
propagators is deiived from Fig. 1 with two derivatives 
w.r.t. the ghost fields and the gluons respectively. The 
propagator flows are coupled sets of integro-differential 
one loop equations for the propagators that also depend 
on vertex functions with up to four legs, see Fig. 2. This 
system of flow equations has been studied in [13, 17, 18, 
19, 20, 21]. In the following we use the numerical solu- 
tions of the FRG equations for the propagators in [13], 
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FIGURE 2. Functional RG equations for the gluon and 
ghost propagator. Filled circles denote dressed propagators and 
empty circles denote dressed vertex functions. Crosses denote 
the cut-off insertion dtR- 

the gluon propagator is displayed in Fig. 3. The vertices 
used for this solution are fully RG-dressed and depend 
on one (symmetric) momentum scale k. The four gluon 
vertex is partially 2PI-resummed and hence includes also 
the sunset diagram with full propagators. If compared 
with similar DSE computations, see [13, 23, 24, 25, 26] 
and references therein, it is in particular the latter prop- 
erty which includes contributions that are contained in 
the two-loop diagrams in the DSEs. These diagrams are 
neglected in most approximations to DSEs (see how- 
ever [27]). This fact is most probably responsible for the 
minor deviations of the DSE propagators from the lat- 
tice results, e.g. [22, 28, 29, 30], in the mid momentum 
regime around the peak of the gluon dressing function 
in Fig. 3, for more details see [13]. In spite of this both 
functional approaches provide results which are in quan- 
titative agreement with the lattice propagators. 
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FIGURE 3. Gluon dressing function: FRG [13], lattice [22] . 



Quark confinement from Yang-Mills propagators 
The Polyakov loop variable L{x), 



— trPexpi^/ <ifAo(f,x), 
A'c Jo 



(3) 



in QCD with A^^ colors and infinitely heavy quarks is re- 
lated to the operator that generates a static quark. The 
trace in (3) is evaluated in the fiindamental representa- 
tion, P stands for path ordering, and )3 = l/T is the in- 
verse temperature. We can interpret the logarithm of the 
expectation value (L) as half of the free energy Fgq of 
a static quark-anti-quark pair at infinite distance. More- 
over, (L) is an order parameter for the center symmetry 
of the gauge group, see e.g. [31]. 

It can be shown that L[{Aq)] also serves as an order pa- 
rameter, [32, 33]: we have L[(Ao)] > {L) with help of the 
Jensen inequality. We also have L[(Ao)] = in the center- 
symmetric phase. Hence it vanishes exactly at T^conf- 

Ac- 
cordingly we have to simply find the solution to the equa- 
tion of motion for Aq: Vyp^[(Ao)] = 0. The effective po- 
tential Vym [^o] is nothing but the effective action Y eval- 
uated for constant fields Aq, Vym[^o] = r[A = Ao,C = 
0, C = 0] . Its flow is governed by the ghost and gluon di- 
agrams in Fig. 1, and hence can be computed solely from 
the (^-dependent) ghost and gluon propagators [32]. In 
other words, VVm being confining for low temperatures 
puts constraints on the behaviour of the ghost and gluon 
propagators in the deep infrared. Loosely speaking, con- 
finement demands a minimal amount of infrared ghost 
dominance [32]. Interestingly, infrared stability in the 
background Landau gauge puts an upper bound on the 
amount of infrared ghost enhancement [34]. Together 
this puts rather tight constraints on the infrared asymp- 
totics of the propagators which are satisfied by the actual 
numerical solutions [34]. Note, however, that the critical 
temperature is insensitive to the deep infrared [32, 35], 
see also [36, 37]. Moreover, confinement is not directly 
sensitive to the size of the coupling. These statements 
hold true in dynamical QCD and extends to other ob- 
servable s [12], 

Inserting the propagators shown in Fig. 3 we get the 
effective potential shown in Fig. 4. The transition tem- 
perature is computed as Tc = 276 ±10 MeV in quantita- 
tive agreement with the lattice result of Tc = 270 MeV, 
e.g. [38]. The order parameter, L[(Ao)] is obtained from 
the minimum of the potential in Fig. 4 and is shown in 
Fig. 5. If compared with lattice results for (L), e.g. [38], 
the Polyakov loop variable L[(Ao)] shows a steeper slope 
and hence a smaller transition region, see also [39], in 
accordance with the Jensen inequality discussed above. 
Note, however, that the results presented in Figs. 4 and 
5 are obtained by neglecting the back-reaction of the Aq- 
fluctuations related to Vym[(Ao)]. These fluctuations re- 
duce the slope of the Polyakov loop variable L[(Ao)]. For 
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FIGURE 4. Polyakov loop effective potential [32]. 

a system at a second order phase transition they carry the 
universal properties of a system and drive the system into 
the symmetric phase. In turn, they are sub-leading for the 
critical temperature of a first order phase transition. This 
is one of the reasons for the good quantitative precision 
of the results for the critical temperature for SU (3) and 
higher SU{N) [35]. 

For SU (2) Yang-MiUs theory we have a second or- 
der phase transition (Ising universality class), and in- 
deed we find Tconf/y/^ = -605 instead of T^-onf/y/a = 
.709 [38], with string tension <7. An alternative compu- 
tation of the order parameters in SU{2) has been done in 
Polyakov gauge. There the flow is completely described 
by the VyMiAo] -fluctuations [33], and the critical temper- 
ature is computed as T^onf/V^ — -69. We also find the 
Ising class critical exponents. More recently, the Landau 
gauge computation has been extended to the VymI^o]- 
fluctuations with Tconf/ a/o' = .705 [40]. A comparison of 
the temperature-dependence of the order parameter in the 
Polyakov gauge and the Landau gauge gives a remark- 
able agreement for all temperatures and provides strong 
support for the gauge independence of our results [33]. 
This concludes the investigation of Yang-Mills theory. 
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FIGURE 5. Polyakov loop [32]. 



QCD, chiral symmetry breaking and dynamical 
hadronisation 

The matter fluctuations due to dynamical quarks are 
encoded in the last two loops in Fig. 1. The quark- 
gluon interaction gives rise to four-fermi terms upon 
integration of a momentum shell with the flow. They 
originate from box diagrams «: d^, leading to 

/ ^v,k{as) [{wwf + {iWJswf] , (4) 

Jx 

where T = (di , (72 , (73 ) . Here we have restricted ourselves 
to the two-flavour case, Nf = 2. Further momentum shell 
integration then also give contributions to the four-fermi 
terms related to (4) as the flow dfX^i^ii^ = dtX^^kiXxi/, oCs). 
Note that for ttj = this resembles a NJL-type flow. The 
four-fermi term in (4) can be conveniently rewritten with 
a Hubbard-Stratonovich transformation as 

\ j 'nl,k{(^^ + + / '^'^ ^^'^^^ ^ + {ifYs^W) ^] , 

(5) 

with m^ J. = h\/{2Xxj,^k) and using the equations of mo- 
tion for (7 and n. Further momentum shell integrations 
generate kinetic terms for the mesonic degrees of free- 
dom a and n and further interaction terms, in particular 
an effective potential Veff.klc^ + which includes the 
mesonic mass in (5). Additionally the four-fermi inter- 
action (4) is re-generated from the quark-gluon interac- 
tion. This already leads to an effective action Tkl'^] with 
^ = (A,C,C, i/A, \^;a,K), and suffices to describe the on- 
set of chiral symmetry breaking, [12, 41, 42, 43, 44]. 
The above structure has been also advocated in [45] on 
the basis of the pure glue results in [33]. An important 
feature of the above bosonisation is that it is not subject 
to double-counting problems. With cut-off terms for the 
mesonic degrees we are led to the flow equation (2) de- 
picted in Fig. 1 . 

In this setting chiral symmetry breaking is monitored 
by the expectation value of a: the Yukawa term in (5) 
can be absorbed into the Dirac term in (1) by 

im^,kWW V {im^,k + hk{(y + iYs^J^)) V- (6) 

Hence, for a non-vanishing expectation value of a this 
term serves as an additional mass term for the quarks. 
We concentrate on the Yukawa term and the effective 
mesonic potential at vanishing pion field, jt = 0, 

fv^s,^[c>^]+ [hk{W)(^- (7) 

Jx Jx 

In the symmetric phase the effective potential Veff^k has 
its minimum at (7 = and the fermionic term vanishes 
on the equation of motion. In the broken phase the min- 
imum is at 5 ^ 0, Vgg[a^] = and the fermionic term 



gives rise to a mass term with mass \hk\o. Indeed, f,i ~ 
\hii\G is an order parameter for the chiral phase transi- 
tion which happens at = 0, that is °°. Note that 
even within this simplified setting this already introduces 
momentum-dependent (non-local) four-fermi couplings. 
Via the couphng to the effective mesonic potential and 
the kinetic terms one also generates (non-local) higher 
order fermionic terms. 

The above approach is systematically improved by 
dynamical hadronisation or rebosonisation |8, 11, 41, 
46]: the re-generated four-fermi interaction can be re- 
absorbed in the Yukawa-interaction in (5). This dynam- 
ically re-adjusts the expansion of the effective action 
in the scale-dependent relevant degrees of freedom and 
guarantees or at least improves the locality of the ex- 
pansion in relevant degrees of freedom. In this way 
the system evolves dynamically from the high tempera- 
ture/large cut-off scale quark-gluon phase to the low tem- 
perature/small cut-off scale hadronic phase, see in partic- 
ular [44]. 

With or without dynamical hadronisation it is a par- 
ticular strength of the present approach that it allows a 
direct access to the physics mechanisms. The role of the 
gauge coupling for chiral symmetry breaking is already 
easily displayed in the setting without (full) dynamical 
hadronisation. In the perturbative regime of QCD the 
four-fermi coupling can be safely put to zero. The afore- 
mentioned box diagrams then generate and successively 
increase the strength of the four-fermi coupling. If it ex- 
ceeds a critical strength, spontaneous chiral symmetry 
breaking is induced similarly to the NJL model. In the 
present QCD approach, this is triggered by the increase 
of the gauge coupling as', for «( < ai crit spontaneous 
chiral symmetry breaking is triggered [12, 41, 42, 43]. 
In turn, if ttj does not exceed the critical couphng a, > 
Cf.v.crit^ the four-fermi coupling never grows big and runs 
into the GauBian infrared fixed point Ay/ = with chiral 
symmetry. We conclude that chiral symmetry breaking is 
primarily driven by the strength of the gauge couphng in 
contradistinction to the deconfinement transition. 



Results: phase structure at vanishing density 

The approach described in the preceding sections is 
now put to work in two flavour QCD in the chiral limit 
and at vanishing density. We map out the confinement- 
deconfinement phase transition with the Polyakov loop 
order parameter L[(Ao)], and the pion decay constant 
f:^ ~ (a). The confinement-deconfinement phase tran- 
sition is also accessible by so-called dual order param- 
eters which are derived from quark correlation functions 
[36, 47, 48, 49, 50] with non-trivial temporal boundary 
conditions for the quark fields in the chosen correlator. In 
[12] it has been shown that such order parameters can be 
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FIGURE 6. Pion decay constant, dual pressure difference 
and Polyakov loop as functions of temperature, Xl = djL, 
XduA = ^7/5 [12]. 



also derived from QCDg at imaginary chemical potential 
/i —IniO/P, for more details see also [51] in these pro- 
ceedings. Here we only remark that imaginary chemical 
potential can be recast as a non-trivial boundary condi- 
tion for the quarks, 



/ 



(8) 

It can be shown that the first Fourier moment in 9 of 
any observable is sensitive to center symmetry [12, 51]. 
Hence, if it does not vanish identically, it is an order pa- 
rameter for the confinement-deconfinement phase tran- 
sition. This allows to define the so-called dual pressure 
(or density) as an order parameter for the confinement- 
deconfinement phase transition (in pure Yang-Mills), 



d0e-2'^"'re[(Ao>e=o,(cT)e=o]. 



(9) 



Here, Tg is the effective action of QCDg, and (Ao)e=o 
and (<7)q=o are the solutions of the equations of motion 
in the physical theory at = 0. The advantage of the 
order parameter (9) is, that its flow dtp is directly related 
to the flow of the effective action, (2). Clearly, this is the 
object least sensitive to a given approximation of F^.. 

The results for the order parameters are displayed in 
Fig. 6 and suggest a relation between the chiral phase 
transition and the broad confinement-deconfinement 
cross-over with a width of approximately 20 MeV, both 
happen at about 180 MeV. The agreement between the 
cross-over temperature derived from the Polyakov loop 
and the dual pressure is remarkable as is the apparent 
close relation of their values for all temperatures. The 
latter has been studied in more detail in [12, 51] and we 
defer the reader to this work. Beside its formal relation 
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FIGURE 7. Pion decay constant as a function of imaginary 
chemical potential and temperature [12]. 

it constitutes a non-trivial consistency check of the 
approximation used as the Polyakov loop is dominated 
by gluonic fluctuations whereas the dual pressure is 
dominated by matter fluctuations. Finally we would 
like to remark that qualitatively our findings compare 
well with the lattice findings for 2h-1 flavours, see e.g. 
[52,53,54, 55,56]. 



Results: phase structure at finite density 

A first step towards non-vanishing density is the inclu- 
sion of an imaginary chemical potential. It allows us to 
compute the dual order parameters such as the dual pres- 
sure (9). We also can collect indirect information about 
real chemical potential by continuation. In specific cases 
this may allow to fix the phase structure at real chemical 
potential [57]. 

The pion decay constant is displayed in Fig. 7. To- 
gether with the Polyakov loop it leads to the phase struc- 
ture displayed in Fig. 8. We conclude that the close 
relation between the chiral phase transition and the 
confinement-deconfinement cross-over persists at imagi- 
nary chemical potential. This is in agreement with related 
lattice computations at relatively heavy quark masses, 
e.g. [4, 58]. In contradistinction, in model computations 
this is only seen within an adjustment of the coupling 
of eight-fermi interaction [59]. In the present QCD ap- 
proach these terms with their QCD-induced coupling are 
generated by the flow, no adjustment is required. This 
pattern allows us to access the interesting question of the 
nature of the Roberge- Weiss endpoint. This is important 
for the extension of the computations at imaginary chem- 
ical potential to real chemical potential, [57, 60], and is 
currently under investigation in the present approach. 
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FIGURE 8. QCD phase diagram at imaginary cliemical po- 
tential. Tlie grey band represents the width of Xl. Black dots in- 
dicate the endpoints of the Polyakov loop RW transitions [12]. 




FIGURE 9. Pure glue part of the Polyakov loop potential in 
two-flavour QCD. 



Finally we consider real chemical potential. As has 
been stressed before, functional approaches do not suffer 
from the sign problem that is so virulent within lattice 
simulations at real chemical potential. A FRG study of 
the chiral phase boundary in one-flavour QCD has been 
put forward in [44] . The results include in particular the 
curvature of the chiral phase boundary which is consis- 
tent with lattice results for small chemical potential. 

First computations within two-flavour QCD at real 
chemical potential have also been performed. Most no- 
tably the two critical temperatures stay close to each 
other; indeed, the confinement-deconfinement tempera- 
ture tends to get smaller in comparison to the chiral one. 
This is not seen in most computations with Polyakov 
loop-extended models and is related to the missing back- 
coupling of the matter sector into the glue sector in these 
models. To see this more clearly we use the fact that the 
above-mentioned models, e.g. [61, 62, 63, 59] and their 
non-local versions, e.g. [64], can be interpreted as spe- 
cific approximations of the fully dynamical QCD flow in 
[12], for more details on this relation see also [65]. We 
first remark that the matter sector including fluctuations 
is given by the last two matter diagrams in Fig. 1 . Switch- 
ing off the glue contributions to the matter propagators 
leads to general quark-meson models which have been 
studied intensively beyond mean field with the FRG, for 
reviews see e.g. [66, 67, 68]. 

In turn, the glue part of QCD is encoded in the first 
two, pure glue, diagrams. Switching off the matter con- 
tributions to the glue dynamics reduces the glue part to 
Yang-Mills theory, and the Polyakov loop potential re- 
duces to that of Yang-Mills theory. In the Polyakov loop- 
extended models, the coupling of these two, now decou- 
pled, sectors is reintroduced via the Aq- (or Polyakov 
loop L) and a-background dependence of the quark loop. 
Note that such an approximation of the QCD flow in 



Fig. 1 still involves a fluctuating matter sector. This de- 
scribes these models beyond mean field as has been stud- 
ied in [65, 69, 70]. 

The above connection of the Polyakov loop-extended 
models to the present QCD approach can be used to qual- 
itatively improve these models: the full QCD effective 
potential (or grand potential) Vqcd is obtained by eval- 
uating the integrated flow in Fig. 1 for constant Aq- and 
a-fields [12], 

Vqcd = Vgiue [Aq , + Vquaik [^0 ■ cr] + Vmeson [^0 1 c] • 

(10) 

The glue potential V^iue is computed from the ghost 
and gluon loop in Fig. 1, Vquark is computed from the 
quark loop, and Vmeson is computed from the mesonic 
loop. Vquark + Vmeson are the direct matter contributions to 
the full Polyakov loop potential. We emphasise that the 
propagators in this computation are that of fully-coupled 
QCD with dynamical quarks and mesons [12]. 

We have already remarked that Vgiue is approximated 
by the Yang-Mills potential VVm in the Polyakov loop- 
extended models. With Vgiue from the QCD flows in [12] 
this approximation can be resolved which will be de- 
tailed elsewhere. Here we discuss how to properly ad- 
just the parameters in the Yang-Mills potential shown 
in Fig. 4 by comparing it to the full glue potential. For 
this purpose we evaluate Vgiue on the solution Ob of the 
equations of motion, Vgiue [Ao] — Vgiue [^o, <7o] as shown 
in Fig. 9. It turns out that the related Polyakov loop com- 
pares well with the Yang-Mills one in Fig. 5 if both are 
plotted as functions of T /Tc. Indeed, the form of the pure 
glue potential is qualitatively unchanged in comparison 
to the Yang-Mills potential, and still shows a first order 
phase transition. The latter has to be taken with a grain 
of salt as we have dropped the center-breaking, though 
sub-leading, Vg^^p-terms in the gluon propagator. In any 
case, the transition temperature is significantly reduced. 
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FIGURE 10. Chiral and deconfinement phase diagram in 
the PQM model. The (grey) band corresponds to the width of 
rf*/rfr at 80% of its peak height [65]. 

from Tym = 276 MeV to Tgiue « 205 MeV. We conclude 
that the full glue potential is modeled well by the Yang- 
Mills potential with a reduced critical temperature. We 
also remark that the above procedure allows in general 
to adjust the sensitive parameters of the Polyakov loop- 
extended models with QCD-input. Most importantly, this 
allows for a systematic improvement of these models to- 
wards full QCD. 

The above analysis confirms quantitatively the phe- 
nomenological HTL/HDL estimate in [61]. Using the lat- 
ter also for larger chemical potential as well as utilising 
FRG flows for the matter fluctuations in the presence of a 
Polyakov loop we are led to the phase diagram displayed 
in Fig. 10, see [65]. At large densities or large chemi- 
cal potential the present approximation without baryons 
certainly is not trustworthy, in particular in the hadronic 
phase. The important inclusion of baryons shall be put 
forward with dynamical hadronisation. Nonetheless the 
combined results from both the model computations and 
the QCD flows constrains the position of the critical 
point; a conservative estimate leads to ii/T > 2/3. 

Conclusions 

In summary we have put forward a quantitative func- 
tional approach to the phase diagram of QCD at finite 
temperature and density. So far fresh insight has been 
gained in Yang-Mills theory, in one-flavour and in two- 
favour QCD at finite temperature and density. The results 
compare well to lattice computations without the fix- 
ing of additional parameters. However, the approach also 
provides results beyond the applicability range of lattice 
simulations, in particular at large density. CuiTently it is 
extended to 2h-1 flavours as well as to baryons. 



Acknowledgments. I thank the organisers for hosting 
such a stimulating conference. I thank J. Braun, A. Eich- 
hom, C. S. Fischer, L. Fister, H. Gies, L. M. Haas, 
T. K. Herbst, A. Maas, F. Marhauser, B.-J. Schaefer, 
F. Spallek and L. von Smekal for discussions and a pleas- 
ant collaboration on the work presented here. This work 
is supported by Helmholtz Alliance HA216/ EMMI. 

REFERENCES 

1. p. Braun-Munzinger, K. Redlich, and J. Stachel (2003), 
nucl-th/0304013. 

2. R. D. Pisarski, and F. Wilczek, Phys. Rev. D29, 338-341 
(1984). 

3. Y. Aoki, Z. Fodor, S. Katz, and K. Szabo, Phys.Lett. 
B643, 46-54 (2006), hep-lat/0609068. 

4. P. de Forcrand, and O. Philipsen, Niicl. Phys. B642, 
290-306 (2002), hep-lat/0205016. 

5. R. Alkofer, et al. (2010), 1012.3192. 

6. C. Wetterich, Phys. Lett. B301, 90-94 (1993). 

7. D. R Litim, and J. M. Pawlowski (1998), 
hep-th/9901053. 

8. J. M. Pawlowski, Annals Phys. 322, 2831-2915 (2007), 
hep-th/0512261. 

9. H. Gies (2006), hep-ph/0 61114 6. 

10. Y. Igarashi, K. Itoh, and H. Sonoda, Prog. Theor. Phys. 
Suppl. 181, 1-166 (2010), 0909.0327. 

11. H. Gies, and C. Wetterich, Phys. Rev. D65, 065001 
(2002), hep-th/0107221. 

12. J. Braun, L. M. Haas, F. Marhauser, and J. M. Pawlowski, 
Phys. Rev Lett. In Press (201 1), 9 8 . 8 . 

13. C. S. Fischer, A. Maas, and J. M. Pawlowski, Annals 
Phys. 324, 2408-2437 (2009), 0810.1987. 

14. L. von Smekal (2008), 0812.0654. 

15. M. Leder, J. M. Pawlowski, H. Reinhardt, and A. Weber 
(2010), 1006.5710. 

16. H. Reinhardt, M. Leder, J. M. Pawlowski, and A. Weber 
(2010), 1011 . 3462. 

17. U. Ellwanger, M. Hirsch, and A. Weber, Z Phys. €69, 
687-698 (1996), hep-th/9506019. 

18. U. Ellwanger, M. Hirsch, and A. Weber, Eur. Phys. J. CI, 
563-578 (1998), hep-ph/9505468. 

19. B. Bergerhoff, and C. Wetterich, Phys. Rev. D57, 
1591-1604(1998), hep-ph/9708425. 

20. J. M. Pawlowski, D. F. Litim, S. Nedelko, and 

L. von Smekal, Phys. Rev Lett 93, 152002 (2004), 
hep-th/0312324. 

21. C. S. Fischer, and H. Gies, J HEP 10, 048 (2004), 
hep-ph/0408089. 

22. A. Sternbeck, E. M. Ilgenfritz, M. MuUer-Preussker, 
A. Schiller, and I. L. Bogolubsky, PoS LAT2006, 076 
(2006), hep- lat/ 610 053. 

23. R. Alkofer, and L. von Smekal, Phys.Rept. 353, 281 
(2001), hep-ph/0007355. 

24. C. S. Fischer, J. Phys. G32, R253-R291 (2006), 
hep-ph/0605173. 

25. D. Binosi, and J. Papavassiliou, Phys. Rept. 479, 1-152 
(2009), 0909.2536. 

26. R. Alkofer, C. S. Fischer, M. Q. Huber, F J. Llanes- 
Estrada, and K. Schwenzer, PoS CONFINEMENTS, 019 
(2008), 0812 .2896. 



27. J. C. Bloch, Few Body Syst. 33, 111-152 (2003), 

hep-ph/0303125. 

28. O. Oliveira, and P. Silva, Eur.Phys.J. C62, 525-534 
(2009), 0705 . 0964. 

29. A. Cucchieri, and T. Mendes, PoS LAT2007, 297 (2007), 
0710 . 0412. 

30. I. Bogolubsky, E. Ilgenfritz, M. Muller-Preussker, and 
A. Stembeck, PoS LAT2007, 290 (2007), 0710.1968. 

31. J. Greensite, Prog.Part.Nucl.Phys. 51, 1 (2003), 

hep-lat/0301023. 

32. J. Braun, H. Gies, and J. M. Pawlowski, Phys. Lett. B684, 

262-267 (2010), 0708.2413. 

33. F. Marhauser, and J. M. Pawlowski (2008), 0812.1144. 

34. A. Eichhorn, H. Gies, and J. M. Pawlowski (2010), 
1010 .2153. 

35. J. Braun, A. Eichhom, H. Gies, and J. M. Pawlowski, Eur. 
Phys. J. C70, 689-702 (2010), 1007.2619. 

36. C. S. Fischer, Phys. Rev. Lett. 103, 052003 (2009), 
0904.2700. 

37. C. S. Fischer, A. Maas, and J. A. Muller, Eur. Phys. J. 
C68, 165-181 (2010), 1003.1960. 

38. B. Lucini, M. Taper, and U. Wenger, JHEP 01, 061 
(2004), hep-lat/0307017. 

39. A. Dumitru, Y. Guo, Y. Hidaka, C. P. K. Altes, and R. D. 
Pisarski (2010), 1011.3820. 

40. J. Braun, H. Gies, J. M. Pawlowski, and F. SpaUek (2010), 
in preparation. 

41. H. Gies, and C. Wetterich, Phys. Rev. D69, 025001 
(2004), hep-th/0209183. 

42. H. Gies, and J. Jaeckel, Eur. Phys. J. C46, 433-438 
(2006), hep-ph/0507171. 

43. J. Braun, and H. Gies, JHEP 06, 024 (2006), 
hep-ph/0602226. 

44. J. Braun, Eur. Phys. J. C64, 459-482 (2009), 
0810.1727. 

45. K.-I. Kondo, Phys. Rev. D82, 065024 (2010), 
1005.0314. 

46. S. Floerchinger, and C. Wetterich, Phys. Lett. B680, 
371-376 (2009), 0905.0915. 

47. C. Gattringer, Phys. Rev. Lett. 97, 032003 (2006), 
hep-lat/060501B. 

48. F. Synatschke. A. Wipf, and C. Wozar, Phys. Rev. D75, 
114003 (2007), hep-lat/0703018. 

49. E. Bilgici, F. Bruckmann, C. Gattringer, and C. Hagen, 
Phys. Rev. D77, 094007 (2008), 0801.4051. 

50. B. Zhang, F. Bruckmann, C. Gattringer, Z. Fodor, and 
K. K. Szabo (2010), 1012.2314. 

51. L. M. Haas, J. Braun, and J. M. Pawlowski (2010), 
1012 .4735. 

52. Y. Aoki, S. Borsanyi, S. Durr, Z. Fodor, S. D. Katz, et al., 
JHEP 0906, 088 (2009), 0903.4155. 

53. M. Cheng, S. Ejiri, P. Hegde, F. Karsch, O. Kaczmarek, 
ct al., Phys.Rev D81, 054504 (2010), 0911 . 2215. 

54. S. Borsanyi, Z. Fodor, C. Hoelbling, S. D. Katz, S. Krieg, 
etal. (2010), 1011 .4230. 

55. A. Bazavov, and R Petreczky (2010), 1 1 2 . 1 2 5 7 . 

56. K.Kanaya (2010), 1012.4235. 

57. P. de Forcrand, and O. Philipsen, Phys. Rev. Lett. 105, 
152001 (2010), 1004.3144. 

58. M. D'Elia, and M.-R Lombardo, Phys. Rev. D67, 014505 
(2003), hep-lat/0209146. 

59. Y. Sakai, K. Kashiwa, H. Kouno, and M. Yahiro, Phys.Rev. 
D77, 051901 (2008), 0801 .0034. 



60. M. D'EUa, and F. Sanfilippo, Phys. Rev. D80, 111501 
(2009), 90 9. 0254. 

61. B.-J. Schaefer, J. M. Pawlowski, and J. Wambach, Phys. 
Rev. D76, 074023 (2007), 0704.3234. 

62. K. Fukushima, Phys.Lett. B591, 277-284 (2004), 
hep-ph/0310121. 

63. C. Ratti, M. A. Thaler, and W. Weise, Phys.Rev. D73, 
014019 (2006), hep-ph/ 050 6234. 

64. T. Hell, S. Roessner, M. Cristoforetti, and W. Weise, 
Phys.Rev D79, 014022 (2009), 0810.1099. 

65. T. K. Herbst, J. M. Pawlowski, and B.-J. Schaefer, Physics 
Letters B In Press (20 1 0), 1 8 . 8 1 . 

66. J. Berges, N. Tetradis, and C. Wetterich, Phys.Rept. 363, 
223-386 (2002), hep-ph/0005122. 

67. B.-J. Schaefer, and J. Wambach, Nucl. Phys. A757, 
479^92 (2005),nucl-th/0403039. 

68. B.-J. Schaefer, and J. Wambach, Phys. Part. Nucl. 39, 
1025-1032 (2008), hep-ph/ 611191. 

69. V. Skokov, B. Stokic, B. Friman, and K. Redlich, 
Phys.Rev C82, 015206 (2010), 1004.2665. 

70. V. Skokov, B. Friman, and K. Redlich (2010), 
1008.4570. 



